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1 Introduction 

An alternative scheme for quantization based on stochastic averages has been presented 
several years ago by Parisi and Wu |l|], see refs. 0, 0] for comprehensive reviews and 
referencing. The main idea of "stochastic quantization" in configuration space - in this 
paper we will exclusively discuss the configuration space quantization - is to view Eu- 
clidean quantum (field) theory as the equilibrium limit of a statistical system coupled 
to a thermal reservoir. This system evolves in a new additional time direction which is 
called stochastic time until it reaches the equilibrium limit for infinite stochastic time. 
The coupling to the heat reservoir is simulated by means of a stochastic noise field - or 
rather by the mathematically well defined Wiener process - which forces the original Eu- 
clidean field to wander randomly. In the equilibrium limit the stochastic averages become 
identical to ordinary Euclidean vacuum expectation values. 

There are two equivalent formulations of stochastic quantization due to the general 
properties of stochastic processes: In one formulation all fields have an additional depen- 
dence on stochastic time. Their stochastic time evolution is determined by a Langevin 
equation which has a drift term constructed from the gradient of the classical action of 
the system. The expectation values of observables are obtained by ensemble averages over 
the Wiener measure. 

Corresponding to this Langevin equation is an equivalent diffusion process which is 
defined in terms of a Fokker Planck equation for the probability distribution character- 
izing the stochastic evolution of the system. Now expectation values of the observables 
are defined by functionally integrating over them with the stochastic time dependent 
Fokker-Planck probability distribution. 

In several non-gauge models (see especially ^ for a mathematically rigorous treatment 
of scalar field theory) one can verify that in the infinite stochastic time limit this Fokker- 
Planck probability distribution converges to the standard Euclidean configuration space 
path integral density. 

Over the past years the stochastic quantization scheme has grown into a useful tool in 
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several areas of quantum field theory and a large number of new results and generalized 
schemes have appeared ^ : In specific we would like to recall the applications to gauge 
fields (see below), fermionic fields, the studies of anomalies, investigations of supersym- 
metric models and the regularization and renormalization program; last but not the least 
the numerical simulations became an important non perturbative application of stochas- 
tic quantization; stochastic techniques could be applied to Monte-Carlo simulations by 
means of the so called hybrid algorithms. 

One might attach great importance especially to the numerical applications of the 
stochastic quantization scheme while supposing that its formal field theoretic develop- 
ments have been completed. It is our conviction, however, that the method of stochastic 
quantization can be distinguished by its fundamental concepts and therefore seems still 
relevant for actual research. In this connection we refer to the quantization of gauge 
theories: 

In configuration space the most prominent quantization scheme is provided by the 
path integral using the Faddeev-Popov procedure 0, which guarantees the unitarity of 
the corresponding Feynman amplitudes. However, due to the Gribov ambiguity |^, 
which is related to the fact that the topological nontrivial structure of the space of gauge 
orbits prevents the existence of a global gauge fixing, the path integral quantization is 
incomplete beyond the perturbation theory. On the perturbative level this does not lead 
to difficulties since one stays within the local region where the gauge fixing is well defined. 
There have been many attempts to deal with these nonperturbative aspects in the path 
integral [|10], |ri|, 0, [1^, but a profound solution has not been found yet. 

Much hope has been put forward to gain new insights for a correct nonperturbative 
path integral formulation of gauge theories from the stochastic quantization point of 
view . The stochastic quantization scheme of gauge theories in configuration space 

has been developed apart of the above dual structure of a Langevin and Fokker-Planck 
formulation in two related ways: 

a) Parisi and Wu [I]] proposed the quantization in terms of Langevin equations with- 
out gauge fixing terms and without Faddeev-Popov ghosts. Since the pure action 
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of the gauge model remains invariant under gauge transformations, the associated 
drift term of the Langevin equation acts orthogonally to the gauge orbits. As a 
consequence unbounded diffusion along the orbits takes place so that potential di- 
vergencies are arising for the expectation values of general observables. It is argued, 
however, that these contributions are canceling out among themselves in expecta- 
tion values of gauge invariant observables. 

This Parisi-Wu stochastic quantization method of gauge theories preserves the full 
gauge symmetry and provides a manifestly gauge covariant kind of quantization. 
Due the absence of gauge fixing it was claimed to be valid also in the nonpertur- 
bative regime without the Gribov ambiguity. 



b) In the approach of Zwanziger |]T^ the drift term of the stochastic process is modified 
by the addition of a gauge fixing force tangent to the gauge orbits. It provides damp- 
ing for the gauge modes' diffusion along the orbits leaving, however, unchanged all 
gauge invariant expectation values. This "stochastic gauge fixing" procedure allows 
to directly resort to a well defined Fokker-Planck formulation so that powerful func- 
tional techniques can be applied. Again no Faddev-Popov ghosts are needed and 
no Gribov ambiguity was claimed to exist neither. 

A fundamental question arises how stochastic quantization - if at all - compares with 
the conventional quantization schemes in the case of gauge theories. A great amount 
of work shows agreement at the perturbative level up to some order in the coupling 
constants, e.g. for Yang-Mills theory up to two loops The question of equivalence, 
however, becomes subtle at the perturbative level to any order and turns into a hot 
subject at the nonperturbative level. 

Unfortunately - due to the unbounded diffusion along the gauge orbits - the Parisi -Wu 
scheme for gauge theories lacks an immediate Fokker-Planck formulation and equilibrium 
properties are hard to extract. This comes from the technical difficulty to solve exactly 
the corresponding Langevin equations; no general equivalence proofs of any nontrivial 
gauge model have been established. 
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On the other hand using Zwanziger's scheme it is necessary to 

a) rigorously control the required stabilizing effect along the gauge orbits 

b) prove convergence to equilibrium of the Fokker-Planck probability distribution 

c) explicitly find out what this limit is. 

With respect to these requirements according to our opinion again no satisfactory 
nonperturbative treatment of any nontrivial gauge model has been given (see however 
the formal procedures |]TU|, PU|). 

In view of the above difficulties it is not surprising that the interest in the stochastic 
quantization scheme - despite of its conceptional beauties - has gone down remarkably 
over the last years . 

It seems first of all helpful to choose a gauge model which is as simple as possible, 
containing, however, relevant physical and mathematical structure. 

From a pure mathematical viewpoint, a gauge model coupled to matter fields can 
be formulated in terms of fiber bundle theory. The fields form an infinite dimensional 
configuration space on which an action of the symmetry group of the model is given. 
Reducing to the true degrees of freedom corresponds to projecting the theory down to the 
true configuration space, which is obtained by taking the quotient of the configuration 
space by the gauge group. Depending on the kind of action of the gauge group this 
quotient may have a smooth or an orbifold structure. 

The problem of the orbifold structure of the physical phase space has been addressed 



by several authors (see e.g. p3| , p2| , pT|). We remind, however, that in this paper we are 
restricting ourselves to just a configuration space discussion. In this respect we would like 
to quote relations with inequivalent quantizations in two dimensional Yang-Mills theory 
on a cylinder by [^, ^ and the analysis of the singularity structure of the Yang-Mills 
configuration space by [^]. 

In order to gain new insights into the stochastic quantization scheme of gauge the- 
ories it seems important to analyze how these geometrical concepts are encoded in the 
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corresponding Langevin or Fokker-Planck equation, respectively. More precisely, we want 
to study the evolution properties of the stochastic process separately along the true con- 
figuration space as well as along the gauge group. In order to avoid orbifold singularities 
and Gribov ambiguities it is appropriate to study such a gauge model in which 

i) the gauge group acts freely, that is without fixed points. 

ii) the separation into gauge independent and gauge dependent degrees of freedom is 
globally possible. 



It is well known (see e.g. [^) that under these circumstances the true configuration 
space is a smooth manifold. Technically, this separation corresponds to the choice of a 
particular trivialization of the principal fiber bundle on which the stochastic process is 
formulated. 

For these reasons we decided to choose the so called Helix model, which describes 
the minimal coupling of an abelian gauge field with three bosonic matter fields in + 
1 dimensions. This model was originally proposed by deWit and was investigated 
intensively within the Hamiltonian framework by Kuchar [E^. Recently the helix model 



came to new life again |13| in the course of studies on problems with gauge fixing and 



the Gribov ambiguity; it was further analyzed within the BRST approach by [jT4| . 



To clarify we remark that although globally the helix model is free from topological 
gauge fixing obstructions one may nevertheless find specific gauge choices which suffer 
from Gribov ambiguities. The helix model therefore serves as an ideal testing ground for 
treatments of "bad" gauge choices and their associated Gribov problems. 

In the present paper (see also ^U, we perform the stochastic quantization pro- 



cedure by exploiting the globally unproblematic structure of the helix model. We will 
come back to the discussion of the above - artificial - complications of the model in a 
separate publication; there we will be analyzing gauge choices where one is forced to work 
on coordinate patches . 



We present in this paper a refined stochastic gauge fixing procedure and prove non- 
perturbatively the equivalence of the stochastically quantized helix model to the known 
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path integral formulation. We propose in fact a modified stochastic process where not 
only the drift term is changed but also the Wiener process itself. These modifications are 
done in such a way that a proper damping of the gauge modes along the orbits is possible, 
that all gauge invariant expectation values remain unchanged and that the equilibrium 
Fokker-Planck distribution can be determined; it is seen to agree with the known path 



integral density of the helix model |^ 



Although our generalized scheme of stochastic gauge fixing is worked out here for the 
helix model we will give arguments for its general applicability. 

We start in section 2) with a review on basic facts of Ito's stochastic calculus and 
recall in section 3) essentials of the Parisi-Wu stochastic quantization scheme of non-gauge 
models. 

A discussion on the geometric structure of the helix model in configuration space is 
given in section 4). 

In section 5) we present an overview of the stochastic quantization of gauge models 
and derive our generahzation of the stochastic gauge fixing procedure. 

The damping of the gauge degrees of freedom is discussed in section 6). 

In section 7) we derive the equilibrium Fokker-Planck distribution of the helix model 
in a nonperturbative manner and show that is agreeing with the known path integral 
density of this model. 

Section 8) deals with geometrical aspects of the generalized stochastic gauge fixing 
scheme. 

We present an outlook in section 9) and add a short appendix for technical details. 

2 A short review of the Ito Stochastic Calculus 



In this section we review basic facts of the Ito stochastic calculus (see e.g. |^3|], 0], ||35| ) 
assuming that the reader is already familiar with the concepts of stochastic processes. 

A set of Langevin equations for m scalar fields s), i = 1, ..,m, which depend on 
a real coordinate t (which in the following sections will be interpreted as Euclidean time) 
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and the stochastic time coordinate s 



d¥{t, s) = a'{t, s)ds + [ dt'h\{t, t'; s)dW''{t', s) 



(2.1) 



is an abreviation for the associated set of stochastic integral equations 



ra'{t,s)da+ f dt' r b\{t,t']a)dW''{t',a) 




(2.2) 



The a* and the b\ are given functions of the fields which are defined at just the same 
stochastic time s. As a consequence the a' and the 6'^ are local in the stochastic time s 
but are allowed to generally be nonlocal in the other variable t. Later in our paper we 
will in fact discuss cases where the a' and the b\ imply differentiations and integrations 
of the with respect to t. 

dW'lt, s) denote the increments of an m-dimensional Wiener process W\t, s), which 
have the correlations 



The second integral in (2.2) is a stochastic integral; it is defined in the mean square 
sense as the limit of an infinite Riemann-Stieltjes sum, which, however, carries a depen- 
dence on the choice of the positions of the intermediate points. We adopt Ito's choice to 
define the integral according to the product rule 



and call the corresponding stochastic differential equation (2.1) an Ito-Langevin equation. 

Under certain conditions (Lippschitz condition, growth condition, specific class of 
initial conditions) it can be proven that the solution of the Ito-Langevin equation (2.1) is 
a unique nonanticipating function. A function g{s) is called non-anticipating if for s < s' 
it is statistically independent of W{s) — W{s'). 

Expectation values of arbitrary functions of the solutions $*(i, s) are obtained by 
ensemble averages over the Wiener measure, using the Gaussian behaviour of the Wiener 
process and applying repeatedly the correlations (2.3). 



{{dW\t, s)dW''{t', s))) = 25'^5{t - t')ds. 



(2.3) 



b\{t,t';s)dW^{t\s)^ lim b\{t,t';s)\W^{t,s + /^s) -W^{t,s)\ 

As— >0 



(2.4) 
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One consequence of Ito's choice is the mean value formula 

{{T fmt,a))dW%a)))^0 (2.5) 

J so 

which follows from the non-anticipating nature of the solutions of the Ito-Langevin equa- 
tion; / is an arbitrary (non-anticipating) function of 

Next we recall the so called Ito formula, which states that the differential of an 
arbitrary function(al) / of is given by 



f dt2j^b\{tuh;s)dW\h,s)\ 



ds 
(2.6) 



Here it is assumed that / is twice partially continuosly differentiablc with respect to 

Equivalently to (2.6) we can formally define multiplication rules for the differentials 
dW' 

dW'{t,s)dW''{t',s) = 25'^5{t-t')ds 

dW\t,sf+^ = AT = 1,2,... (2.7) 

dW\t,s)ds = 0, 

which - more precisely - are given in terms of associated integrals over arbitrary non- 
anticipating functions. From (2.6) or (2.7) amounts, unfortunately, the cumbersome con- 
clusion that in calculating expressions involving infinitesimals ds of the stochastic time 
one must keep all expansions up to the second order in the dW''. 

What usually is meant by Ito's stochastic calculus can be summarized by equations 
(2.4) - (2.7). 

It is of fundamental significance that with the help of Ito's as well as the mean value 
formula one can derive from the Ito-Langevin equation (2.1) the Fokker-Planck equation 
for a (conditional) probability density p 

^*^ = L[$]p[$,.]. (2.8) 
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Here the so called Fokker-Planck operator L is defined in terms of the a' and the ^ 

It should be stressed that in a Fokker-Planck formulation the fields - and consequently 
also the a* and the b^k - are stochastic time independent; all the stochastic time depen- 
dence is coming from the probability density p itself. We therefore omit to denote a 
s-dependence of as well as of a* and 6*^. 

Due to the fundamental relationship between the Ito-Langevin (2.1) and the Fokker- 
Planck (2.8, 2.9) scheme the expectation values of arbitrary (equal stochastic time) func- 
tions / of the solutions $*(i, s) of the Ito-Langevin equation are coinciding with expec- 
tation values defined correspondingly in terms of the Fokker-Planck density p 

((/($(•, s)))) = {fmm = Jd^ /($(-))p[$, s]. (2.10) 

We remark again that in (2.10) there are involved stochastic time dependent as well as 
stochastic time independent fields, respectively. 

The Fokker-Planck operator L describes the stochastic time evolution of the Fokker- 
Planck probability density; with its help we formally write 

p[$,s] = e^^[*V[*,s = 0]. (2.11) 

In case that we are concerned with the stochastic time evolution of expectation values, 
it is advantageous, however, to introduce in addition to L also its adjoint 

I m{-M^,s] = / L'$(e^^[*lp[$,« = 0])/($(-)) ^Jd^ = 0]e^^'[*l/($(-)). 

(2.12) 

It exphcitly reads in terms of the a' and the b^k 



Lt[$] = ^ dt, a'ih) + J^^ dt2dtsb'j{tu ts)S^^b\{t2, t 



6 

3) 



^ (2.13) 



At the end of this section we remark that we discussed Ito's stochastic calculus for fields 
depending on just the coordinate(s) t (and s), having in mind applications to the helix 
model. It would otherwise have been straightforward to generalize this section to a higher 
dimensional case as well. 
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3 A short review of the Parisi-Wu stochastic quan- 
tization scheme 

In this section we review some of the important features of the stochastic quantization 
scheme when performing the quantization of non-gauge models. We include this section 
for the benfit of the reader to get a self-contained presentation of our paper in order to be 
able to follow closely our subsequent discussions of the stochastic quantization procedure 
in the case of gauge theories. 

We consider a model of m (self-interacting) scalar fields with the positive Eu- 
clidean action denoted by S. We assume that the model has a finite normalization integral 

The basic idea of the stochastic quantization scheme is to interpret the Euclidean path 
integral measure / J as the equilibrium limit of a Fokker-Planck probability 

distribution of a specific stochastic process. 

This stochastic process is defined by an Ito-Langevin equation (2.1) of the following 
form: The drift term a* is taken to be the variation of the action S with respect to the 
fields whereas the matrix 6*^ is chosen most simply as an m-dimensional unit matrix. 

6S 



ds + dW\t,s). (3.1) 

$(t)=<I>(j,s) 



It then follows from (2.9) that the associated Fokker-Planck operator is just 



6S 6 



(3.2) 



We nonperturbatively deduce that the equilibrium limit of the Fokker-Planck distribution 
is formally given by e^^ / J . In this paper we do not discuss the complications 

due to the issue of renormalization, see, however ^ or pn|. 



As a consequence the equal stochastic time expectation values of arbitrary observables 
(calculated either within the Langevin or Fokker-Planck approach, (2.10)) relax in the 
infinite stochastic time limit to the corresponding quantum Green functions 

lim ((/($(-, .)))) = lim (/($(■)))(.) = (/($(■)))■ (3.3) 
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For later application we study now a bijective change of variables of the original fields 
<l>*(t) to new fields \E''^(i(:) With respect to this variable change we introduce 

vielbeins E and their inverses e 



as well as the induced inverse metric G^'^ and the determinant G 



(3.4) 



R 



G = detGf,^{ti,t2). 



(3.5) 



The vielbeins, the inverse metric G'^'^ and the determinant G are generally field dependent 
and are used either in the Fokker-Planck or in the Langevin equations; in the latter case 
we have to explicitly insert the additional s-dependence of the fields. 

With the Ito formula (2.6) the transformed Langevin equations for the new variables 

read 



R 



dt'l E^,(tX) 



5S 



-ds + dW'{t\s) 



which we can recast into 
d<i'^'{t,s) = I dt' 



-G'''{t,t')^J^+ ^ 



(3.6) 



+ E^',{t,t')dW\t\s)] 



^(G-(M')v/G) 



ds 



(3.7) 



*{-)=*(-,^) 

In the last equation we made use of the fact that the covariant divergence of the vielbein 
E is vanishing (we considered the transformations from flat field variables to curved ones) 
so that 



dt 



R 



6E>\(t^ I 6VG 
5^^(t) VG '^5^^(t) 







and 



R 5<^'{t') 



dt'^ — 

R 5^ 



^(G^-(t,OVG) 



(3.8) 



(3.9) 



We interpret (3.7) as a stochastic process taking place in the manifold parametrized by 
the metric G^^, see e.g. we observe the presence of the so called "Ito term", which 
appears as a simple consequence of the Ito formula. 
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We are interested in the transformed form of the Fokker-Planck equation, too. We 
have just to use (2.9) in the case of the new variables starting from the transformed 
Langevin equation (3.7). We immediately obtain the Fokker-Planck operator in the new 
variables as 



m 



dtdt' 



6 



R2 



6S 



dVG 



+ 



G S'^^it') (5^'^(f 



(3.10) 



The formal equilibrium distribution now reads 



Ge- 



(3.11) 



and we observe the appearance of the familiar Jacobian factor corresponding to the 
variable transformation to \E''^. 

Finally we comment on an equivalence relation of stochastic processes in the infinite 
stochastic time limit, which is usually called fluctuation-dissipation theorem, or some- 
times kernel independence of the stochastic process (see e.g. |Q and |3^). More precisely 
we consider Ito-Langevin equations of the form 

6S 5 



d-^^it, s) 



dt' 



-G^''{t,t') 



+ 



+ Ef'i{t,t')dW\t',s)} 



-G'^'M 



ds 



^(■)=^(-,s) 



(3.12) 



which are distinguished from (3.7) with respect to a different form of the Ito term. 
We remark that G^''' is assumed to be positive as above with the determinant G being 
arbitrary. It follows from (3.12) that G'^'^ factorizes in the corresponding Fokker-Planck 
operator 



m 



R2 5^/^(t) 



G'"'it,t') 



6S 



+ 



(3.13) 



so that the equilibrium distribution is given independently of the explicit form of G^'^ by 



equil. 



(3.14) 



13 



4 The geometrical structure of the hehx model 

The helix model is defined by the Lagrange density 

(4.1) 

where the dot denotes time derivation and the fields 

and A{t) are regarded as elements of the function spaces £ — C°°(R, R^) and A — 
C°°(R, R), respectively 

Let Q = C°°(R, R) denote the abelian group of local gauge transformations where 
{9i'92){t) — 9i{'t) + 92{t) gives the group structure. Consider the following transformation 
7?. of ^ on the configuration space £ x A 

n{^,A,g){t) = {R{g{t)Mt),^\t) - g{t),A{t) - g{t)), (4-2) 

where g E Q and 

[ cos^(t) -sin^(t) 

R{9{t)) = . , , 

\^ sm g{t) cos g{t) 

The Lagrange density (4.1) is easily verified to be invariant under these transformations 
and therefore it descends to a function on the quotient space £ Xg A:= {£ xA)/G which 
is obtained by factoring out the gauge degrees of freedom. The generator of infinitesimal 
gauge transformations is given by the vector field Z^{^) with components Z^{t',t) 



dt 



Scp^it) ' ^ ' V(^) V(^) 'SA{t)_ 

- LdtdtW)Z\t',t)Jr^^ (4.3) 

on £^ X ^, where ^ e Lie Q — C°°(R, R) is an arbitrary element of the Lie algebra of the 
gauge group and where we defined 

m = im^m) = {^\t),^'{t),^'{t),A{t)). 
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The transformation properties (4.2) suggest to view the fields A{t) and ip(t) as gauge field 
and matter fields, respectively, so that the helix model can be interpreted as an example 
for a + 1 dimensional gauge theory on the real line R coupled with three matter fields. 
Prom a more mathematical point of view, the space A can be identified with the affine 
space of connections on the trivial principal fiber bundle over the base R with the additive 
group of real numbers as corresponding structure group and the matter fields can then be 
realized as sections of a particular vector bundle with fiber R^ associated to that given 
principal bundle. 

The true configuration space Exg ^ is a smooth manifold which follows from the fact 
that the gauge group Q acts freely on £^ x ^, which means that the equation TZ{0, A, g) — 
{(f, A) for arbitrary fields {^,A) e 6 xA only admits the trivial solution g = 0. Moreover 
the projection n : £ x A ^ S Xg A, with 



where [0, A] denotes the equivalence class of {0, A) with respect to the action TZ, provides 
the structure of a principal fiber bundle with total space £ x A, base manifold £ Xg A 
and structure group Q. 

In this paper we aim at a nonpcrturbative application of the stochastic quantization 
scheme to the helix model. In order to proceed it will turn out to be of great importance 
to rewrite the stochastic process in terms of gauge invariant and gauge dependent fields. 
However, from a mathematical viewpoint, such a rewriting really can be performed if 
and only if the principal fiber bundle n : £ x A ^ £ Xg A is trivializable. In this case 
every field configuration [0, A) can be represented in terms of coordinates parametrizing 
£ Xg A and Q, respectively. 

In order to prove the triviality of the bundle we necessarily have to find a global 
section, that means a map a : £ Xg A ^ £ x A satisfying ir o a — id. We will study two 
specific examples for such global sections. 

First we can easily verify that ai : £ Xg A ^ £ x A given by 



n{0,A)^[0,A], 



(4.4) 



a,{[0,A]) = {R{^')^,O,A-0') 



(4.5) 
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is indeed a global section, corresponding in physical terms to the axial gauge (p^ — 0. 

However wc shall display another global section of the bundle which implies the tem- 
poral gauge ^4 = 0. For this case let us take an arbitrary but fixed function h e C°°(R, R) 
with the property that Jj^ dth{t) — 1 and /^^ dth{tY — 1. The map 02'- S Xg S x A 

given by 

a2iW,A]) = {R{F{^,A))^,^'-F{^,A),0) (4.6) 
where F : S x A^ Q is defined by 



F((f,Am := / dt'h(t')\ t dt"A(t") + (p%t') 
Jr Ut' 



(4.7) 



yields a global section oi n : £ x A ^ £ Xg ^. To see this, we first note that (J2 is well 
defined which follows from F{R{ip, A, g)) = F{0,A) — g. To prove finally that (T2 is a 
section we have to use the fact that F{(f, A) — A. 

For a given section (jj we have a diffeomorphism Xi '■ ^6 A) xQ ^ £ x A given by 

XmA],g)^niamA]),9) (4-8) 

where i = 1,2. Let (^(o-^) : £ x A ^ Q he the function which is uniquely determined by 
the equation 

{<f,A)^n{ai{7r{0,A)),u;^,,){^,A)). (4.9) 

Explicitely we find 

uj^,^){0,A) = a^{[0,A])-p\ (4.10) 

where af{[0,A]) denotes the third component of the i-th section (jj. For the two choices 
(4.5) and (4.6) for the sections (Tj we have 

u;^,,0,A)^-<p' (4.11) 

and 

u;^,,){0,A)^-F{0,A), (4.12) 
respectively. We remark that under a gauge transformation LJ(^ai) transforms simply like 

u^,,){n{0,A,g))=u;^,^){ip,A)+g (4.13) 
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and for later use also mention the important property 

a;(.,)(<7,([(^,A])) = 0. (4.14) 

Now it can be verified straightforwardly that the inverse mapping xT^ : £ x A ^ 
{£ Xg ^) X ^ is given by 

Xl^{0:A)^{[0,Aiu;^^^0,A)) (4.15) 

and it is therefore globally possible to rewrite every field configuration in the appropriate 
manner. In physical terms the triviality of the principal fiber bundle tt : £ x A ^ S Xg A 
means that the helix model does not suffer from a topological Gribov problem. 

However, for further use it is inconvenient to rewrite all fields in terms of equivalence 
classes on S Xg A. Therefore it seems to be more appropriate to introduce coordinates 
which parametrize a space diffeomorphic to the orbit space S Xg A which also may be 
embedded into the original configuration space £ x A. A natural choice is provided by 
the image of the global section a. 

Let S(o..) = im cXi denote the gauge fixing surface in £^ x .4. for a given section cTj 
then we can consider the isomorphic principal fiber bundle 7r(o..) : S x A ^ I^(ai) with 
projection 71(0-^) = (7,077, 

n^,,){0,A)^ai{[0,A]), (4.16) 

total space S x A, and base manifold which can be treated as a submanifold in 

8 X A (see Fig. 1 for these geometrical structures). 

At this point we want to remark that due to (4.14) the fields {0,A) on the gauge 
fixing surface S(o-.) can equivalently be characterized by the implicit condition 

u;^,,){0,A)^O. (4.17) 

For our first choice (4.5) of the global section ai we simply obtain from (4.11), (4.17) 
that (p^ — 0, whereas for (T2 defined in (4.6) we derive 

A^O as well as / dth(tW(t) = (4.18) 
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for h with the properties as given above. 

Corresponding to the bundle ir(^^.^ : S x A ^ ^(o-i) ^ section can be chosen to be the 
identity mapping so that in analogy to (4.8) the bundle trivialization X{ai) '■ ^{ai) x Q ^ 
£ X Ais given by 

Xiad'^mA]),g)^n{amA]):9)- (4-19) 

Corresponding to (4.10) and (4.15) the inverse mapping x^^i) '■ E x A ^ E^g..) x Q 
reads 

Xi.\)(^,A) = (amA]),u;^.MA))- (4-20) 
In subsequent chapters we will use exactly this diffeomorphism S x A = ^{cn) x Q for 
studying the stochastic process associated with the helix model. This result provides the 
way to separate the original process in the configuration space £^x^ into two subprocesses, 
one of which is formulated just in terms of fields {(p, A) e '^{ai) whereas the other process 
is involving gauge dependent fields as well. Since this diffeomorphism explicitly depends 
on the chosen section (T(i), the corresponding processes will be related with the actual 
gauge fixing. 

Let us introduce general coordinates '^{ai) on '^{cn) 

= (s^(.,), *?.,)) = (*U)' "^U^ "^U) 

and a coordinate ^(g. ) parametrizing which we combine into 

Prom (4.20) we trivially deduce 3-dimensional coordinates on E(ct^) by observing that for 
each of the two choices (4.5), (4.6) the 4-vector (Ti([<^, A]) has one vanishing component; 
'^\fJ^) is taken to be a;((^.)(<^, ^4). 

For the particular section ui in (4.5), (4.20) yields 

^ A-^^ (4.21) 
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with inverse 



= (4.22) 

For (72 as defined in (4.6) the corresponding coordinate transformation (4.20) reads 

^(,,) = R{F{^,A))^ 

= ^'-F(0,A) (4.23) 

,4 



with inverse 



= (4-24) 



where 



*^.) W = £ ['^(^ - - M*)M^')] *^.)(^')- (4.25) 
Notice that corresponding to (4.17) and (4.18) the coordinates ^(0-2) on S(cr2) fulfiU 
/r '^^^(^)^((T2)(^) ~ which explains the presence of the projection operator S{t — t') — 
h{t)h{t') acting on ^^^^^ in (4.24) and (4.25). This operator is not needed in (4.23) as 
^\a2) defined in terms of {(p, A) automatically fulfills the required relation . 

In the new coordinates, gauge transformations are given purely as translations, i.e. 
(^(o-j), ^^^^-j) (^(o-.), + g) where g E Q. Correspondingly, the generator for an 
infinitesimal gauge transformation ^ e Lie Q 

= J^/^ (4.26) 

is just the invariant vector field along the gauge group. 
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Let us discuss our findings of tliis cliapter in more detail: We liave seen tliat tlie 
principal fiber bundle associated with the helix model is trivializable and so does not 
suffer from a Gribov problem. To clarify this point we recall that in strict mathematical 
terms (as analyzed by Singer ^ in the case of Yang-Mills theory) one speaks about a 
Gribov problem if the projection from the original configuration space onto the space of 
gauge orbits gives a non-trivializable principal fiber bundle. In this case it is not possible 
to find a separation of any field configuration in terms of gauge invariant and gauge 
variant fields in a global way. Such a separation can be done only locally and as a result 
there may occur ambiguities in the choice of appropriate coordinates on the local patches 
of the orbit space which have nonvoid intersection. 

Our main idea was to select a gauge fixing surface in the original configuration space. 
This choice was naturally given by taking a particular section in the corresponding prin- 
cipal fiber bundle. The gauge fixing surface was then characterized as a set of solution 
of a certain function (4.11) on the configuration space depending, in fact, on the chosen 
section. 



We recall at this place that recently several papers have been published pTS] , |14 
which take the helix model as an (explicitly soluble) example for analyzing the Gribov 
problem. The apparent discrepancy to our result can be traced back to different concepts 
of definitions: There exists of course the possibility to define an arbitrary function on the 
configuration space, consider the set of points where it vanishes and ask if this set is a 
diffeomorphic copy of the gauge orbit space. It is precisely this route of analysis which 
has been followed in the above mentioned papers. Depending on the actual form of the 
given function, it may be possible to find at least a certain subset of solutions which is 
diffeomorphic to a subset of the orbit space. In that case the theory was said to possess 
a Gribov ambiguity. We prefer to interpret the kind of Gribov ambiguity which is talked 
about in this context as being of an artificial kind. 



In a following paper |^ we shall present a stochastic quantization analysis of this 
artificial Gribov problem associated with sections of the bundle TT:SxA^Sxg A 
which can be defined only locally. 
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5 Generalized stochastic gauge fixing 

In this section we review the stochastic quantization procedure of gauge (field) theories 



and present our generahsation of the stochastic gauge fixing procedure of Zwanziger [17 

In order to work out the new features of our modified approach we review at the begin- 
ning the Parisi-Wu stochastic quantization scheme ^ for gauge theories and Zwanziger's 
original gauge fixing procedure. This summary serves as a complement to our general 
introductory remarks in sect. 1) and to our short review in sect. 3) of the stochastic 
quantization of non-gauge models. 

We collectively denote by the pure gauge, as well as matter fields of the given gauge 
model. For reasons of simplicity we choose a notation in correspondance with the helix 
model; in specific we introduce 4 fields depending on the (Euclidean) time coordinate 
t, as the helix model is defined in terms of 4 fields in + 1 dimensions. We remark 
that a generalization to more space time dimension and to more general gauge models 
is immediate. According to the setting of stochastic quantization within the Langevin 
picture the fields are defined to depend additionally on the stochastic time coordinate s 
as well; explicitly this amounts to substitute ^''(t) by $^(s,i(:). 

The crucial point of the Parisi-Wu approach for the gauge theory case is to demand 
that the stochastic time evolution of the fields is given again by a Langevin equation of 
the form 



ds + dW' (5.1) 

<i>{t)=<l>{t,s) 



where S denotes the original (Euclidean) action of the gauge model; it is the unmodified 
bare action without gauge symmetry breaking terms and without accompanying ghost 
field terms; in our case of the helix model it is given by the time integral of the Wick 
rotated form of the Lagrangian (4.1). 

We have ahead remarked in the introductiory section that in the case of a gauge 
theory the stochastic process defined by (5.1) undergoes undamped diffusion and does 
not approach an equilibrium distribution. Related to this fact is that a Fokker-Planck 
formulation for the is not possible because the gauge invariance of the action leads to 
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divergencies in the normalization condition of the Fokker-Planck density. 

We now discuss Zwanziger's modified formulation [|l^ of the Parisi-Wu scheme: The 



stochastic gauge fixing procedure consists in adding an additional drift force to the 
Langevin equation (5.1) which acts tangentially to the gauge orbits. This additional 
term generally can be expressed by the gauge generator and an arbitrary function a; in 
the specific case of the helix model the components of the gauge generator are defined 
by (4.3) so that the modified Langevin equation reads as follows 



d¥{t,s) 



/ dt'Z\t,t')a{t') 



ds + dW' (5.2) 



5^^ it) Jr 

Let us note two remarkable consequences of the extra term in (5.2) 

a) The expectation values of gauge invariant observables remain unchanged for any 
choice of the function a (see below for the explicit demonstration contained in the 
discussion of our generalized stochastic gauge fixing procedure) 

b) For specific choices of the - in principle - arbitrary function a the gauge modes' 
diffusion is damped along the gauge orbits. As a consequence the Fokker-Planck 
density can be normalized and arbitrary expectation values of the fields are relaxing 
exponentially in the stochastic time s to equilibrium values. We remind that this 
situation is in contrast to the Parisi- Wu approach, where for expectation values of 
gauge variant observables no equilibrium values could be attained. 

We present now our generalization of Zwanziger's stochastic gauge fixing procedure 
by adding a specific drift term which not only has tangential components along the gauge 
orbits; in addition we modify the Wiener process itself. Due to this generalization we are 
able to introduce more than just only one extra function a, in fact we additionally intro- 
duce 4 arbitrary functions (3k- Our generalization is done in such a way that expectation 
values of gauge invariant observables again remain untouched. 

Our generalized Langevin equation reads 



d^'{t,s) 
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ds 



+ / dt2 



S'k5{t-t2)+ / dtiZ\t,ti)(3k{tiM) 



R 



*(.)=*(-,s) 



(iW^'(t2,s) 



(5.3) 



here we introduced C'^(^i)^2) as a shorthand notation of 



C\hM) = 25^%(ti,t2) + / dhdUZ\t2,h)l3,{hM)5'"'l3m{tiM)- (5.4) 

We see that the new drift term clearly is not acting tangential to the gauge orbit; its rather 
complicated structure is necessary for leaving unchanged gauge invariant expectation 
values. We prove that this is indeed so. 

Suppose that F is a gauge invariant functional of the fields so that 

5F 



I 



dt 



0. 



(5.5) 



We consider now the difference of the expectation values of F when we first use the 
generalized Langevin equation (5.3) and second the original Parisi-Wu equation (5.1), 
which is just the a = 0, = contribution of (5.3). With the Ito rule (2.7) we find 
straightforwardly 

5F 



{{dF)) - {{dF))l^,^^^^,= {{J^Jt,dt2 



Z'(h,t2)a{t2) 



+ 



L 



dU 



6F 



-Z\h,t2)]c'{t2,h) 



))ds 



(5.6) 



which is zero, as each term vanishes separately due the gauge invariance of F. 

Our next task is to discuss the issue of stochastic gauge fixing in terms of the Fokker- 
Planck equation. Zwanziger originally worked in this scheme to derive his stochastic gauge 
fixing procedure; wc will now generalize his formulation. 

A remarkable advantage of the Fokker-Planck scheme in comparison to the Langevin 
formulation is that equilibrium properties of the system may directly be studied in a 
nonperturbative way. By this we mean that one neither needs to explicitly find the so- 
lution of the stochastic differential equation of the fields nor is it necessary to solve 
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the Fokker-Planck equation for the probabihty distribution p for all stochastic times s. 
What, however, one may directly be able to achieve is to extract from the Fokker-Planck 
equation just the equilibrium distribution. Surely one must be cautious to check that all 
conditions for the stochastic process really to relax to this equilibrium distribution are 



We will demonstrate below that our generalized stochastic gauge fixing formulation 
will indeed allow us to extract in a nonperturbative manner an equilibrium Fokker-Planck 
probability distribution which defines a standard field theory path integral density for 
the helix model. 

First we want to repeat within the Fokker-Planck approach that gauge invariant 
expectation values remain unchanged when our generalized stochastic gauge fixing pro- 
cedure is applied. We have to study the stochastic time development of gauge invariant 
expectation values and therefore are interested in the form of the adjoint Fokker-Planck 
operator . Corresponding to the Langevin equation (5.3) and following the general 
discussion (2.1) and (2.12) we obtain 



We recognize that each of the last two terms is a differential operator acting along the 
gauge orbits so that it vanishes when it is applied to gauge invariant observables. Thus we 
have demonstrated within the Fokker-Planck formulation as well that expectation values 
of gauge invariant observables are not affected by our modifications. 

Comparing with Zwanziger's original stochastic gauge fixing procedure our general- 
ized scheme amounts to adding to the (adjoint) Fokker-Planck operator not only a gauge 
symmetry preserving part which is linear in the field derivatives, but also a gauge sym- 
metry preserving part which contains quadratic derivatives. In other words not only the 
drift term but also the diffusion term of the (adjoint) Fokker-Planck operator gets modi- 



fulfilled. 




+ 



(5.7) 
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fied; the modification of the latter corresponds precisely to the modifcation of the Wiener 
process of the Langevin equation (5.3). 

It is of crucial importance in Zwanziger's approach to draw advantage of the fact 
that the expectation values of gauge invariant quantities are independent of the new 
drift term, which in specific means that they are independent of the function a. One 
may chose now such a specific value for it that the functional integrations along the 
gauge orbits receive proper damping. Equivalently the gauge modes' contributions to 
expectation values are seen to converge rapidly, i.e. exponentially, in stochastic time. 
Only with such a proper choice the Fokker-Planck density can be normalized to unity 
and allows for a true probabilistic interpretation of the stochastic scheme. 

We want to discuss the issue of proper damping of the gauge modes within our gen- 
eralized scheme in the next section. 

6 Damping of the gauge degrees of freedom 

So far we explained our stochastic gauge fixing method in quite general terms, but we 
now want to consider gauge models where the transition to gauge invariant variables at 
the one hand and gauge transforming variables on the other hand globally is possible. 
Geometrically this corresponds to a situation where the configuration space admits a 
structure of a trivialisable principal fiber bundle. The transition to these new variables 
is not unique and corresponds to a specific choice of gauge. 

In the following we specialize further and discuss the helix model in detail. In this way 
we can make use of our geometrical analysis in chapter 4) and concentrate on the two sets 
(4.15, 4.17) of field variables ^(o-^), i = 1, 2, which are related to the gauge fixing surfaces 

= and ^4 = 0, respectively. We remind that in both cases the gauge dependent 
variables are changing under gauge transformations just by a simple shift. 

We want to investigate how the Langevin equations (5.3) are changed when we trans- 
form the original field variables $ into one of the sets of new coordinates ^(o-^). With 
respect to these variable changes we introduce the vielbeins E and their inverses e, as 
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well as the induced inverse metric G'^'^; see (3.4-3.5), as well as the Appendix for expUcit 

expressions. 

With the use of the Ito formula (2.6) the Langevin equations (5.3) are transformed 
in a similar way as (3.7) into 

5S S^^^it) 



(i^^(t,s) 



R 



+ 



ds 



+ f dt2E^',{t,ti)Z\tiM)oi{t: 
Jn 

+ [ dt2E\{t,h) \s\S{tl-t2)+ f dtsZ%,ts)(3k{t3,t2) 



dW''{t2,s) 



R3 



dt2dtsdt4 



-{E^i{t,U)Z\U,ti)) 



ati,t2) 



(6.1) 



*(.)=*(-,s) 



5^''{t2) 

Fortunately the last term in this somewhat involved looking expression vanishes. This 
follows from 

Z\t,t') ^ e\{t,t') (6.2) 

and the fact that e is the inverse to E. The expression (6.2) can be derived by consid- 
ering that the new coodinates ^ are gauge invariant except for which, however, is 
transforming only by a linear shift, see the discussion below (4.13). Therefore under an 
infinitesimal gauge transformation ^ e Lie Q the fields '^^{t) transform according to the 
Lie-derivation along the vector field Z^, 

z^{^^{t)) = l^^dhdt2mz%M^^^ 



/R2 

so that (6.2) can be obtained. We also have 



dhdt2^{tl)Z\h,t2)E>'i{t,ti) 



L 



(6.3) 



(6.4) 



which is satisfied by both of the above sets of coordinates ^^^r) (see the Appendix) so 
that the Ito term takes the special form 



I dh 



[ dhdt2E\{t2,h)- 
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6 



Thus we finally find the Langevin equations for the new variables \Ef(o-) 



(6.5) 



d'^>'{t,s) 



R 



dti 



6S 



+ [E^,{t,t,) + 6%(3,{t,h)]dW\h,s)} 



+ (5^(5(t-ti)a(ti) 



ds 



(6.6) 



We recognize that the terms due to the generalized gauge fixing procedure - the terms 
containing the arbitrary functions a and jSk - are not affecting the Langevin equations of 
the gauge independent fields ^(o-); /i = 1, 2, 3 

6S 6G^''{ti,ty 



c/^^(t,s) 



R 



dti 



+ EMt,t,)dW\h,s)} 



ds 



(6.7) 



*{-)=*{-,^) 

It is important to remark that these Langevin equations are closed among themselves 
and can be solved independently of the gauge dependent variable \['^. This may look 
surprising, as the vielbein E generally is a gauge dependent matrix {G^'^, however, is 
gauge independent by direct inspection). The claimed gauge independence of the above 
Langevin equations follows because 

i) the vielbein E can be split into the product of a gauge invariant matrix and a pure 
rotation matrix which contains all the gauge dependence (see (A. 3) and (A. 7)) 



ii) the Ito-Langevin equations are invariant under rotations of the vielbeins (see 
and references therein), so that the gauge dependent part of E can just be rotated 
away. 

iii) The gauge independence of G^'^ follows generally from (3.5) as a direct consequence 
of the product form of E. 

At this point it seems appropriate to notice that equation (6.7) can also be obtained 
within the pure Parisi-Wu formulation of gauge theories, i.e. without introducing any 
stochastic gauge fixing terms. In this case the original Langevin equation (5.1) constitutes 
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the starting point for introducing the new coordinates all the steps concerning the 
application of the Ito formula are the same as above. The formulation (6.7) in terms of 
only gauge invariant fields is precisely in the spirit of the original paper of Parisi-Wu [|^ 
where the transition to gauge invariant variables and the discarding of the gauge variant 
variables was suggested. In our case we are, however, particulary interested in connecting 
the stochastic formalism to the path integral scheme. We therefore prefer to consider the 
stochastic time evolution of all the variables of the gauge model and continue within the 
generalized stochastic gauge fixing approach. 

In order to derive the Langevin equation of the gauge dependent field we use 



This follows from (6.4) and the fact that \I'^ is linear in see (4.15) and (4.17), as well 
as (4.6). We obtain 



We note that on the right hand side of (6.9) no gauge dependent terms are present apart 
of the possibility that the so far arbitrary functions a and are changing under gauge 
transformations. Let us suppose, however, for a moment that a is gauge invariant. In 
this case the drift term of (6.9) is gauge invariant as well and remains unchanged for 
a stochastic evolution along the gauge orbits, i.e. along the gauge transforming variable 

keeping the other variables fixed. For this stochastic process the drift term of the 
Langevin equation of appears like a constant; it gives rise to a diffusion behaviour 
and will never allow relaxation to an equilibrium distribution. 

But we can convert this diffusion into a well damped stochastic process by conveniently 
choosing a gauge dependent a: We recall that due to the general philosophy of stochastic 
gauge fixing, which we have formulated in the previous section, every choice of a (as well 
as of /3fc) is permissible. 

Out of the possible choices the simplest Langevin equation is obtained if a subtracts 
off the above gauge invariant drift term and adds a term proportional to ^^'^ itself, the 






(6.9) 
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proportionality factor being built from a gauge invariant and strictly positive integral 
kernel 7 



a{t) = j^dt' [G'^M^^^^ - l{t,m\t') 



(6.10) 



We can suggestively rewrite this damping term by introducing the total action jS'tot 



2 Jn 



(6.11) 



so that the Langevin equation becomes 



+ 



(6.12) 



Let us define a matrix G by 



G^"(ti,t2) = G'^'(ti,t2) fi,u = 1,2,3 



(6.13) 



Notice that at this moment it is not obvious that G as defined above is a positive matrix 
and therefore satisfies the properties of a metric; we will prove it, however, in the next 
section for a special choice of 13k and 7. Since from (6.8) 

the Langevin equations (6.6) can be rewritten in terms of the total action ^tot and the 
new metric G as follows 

5S,,, , 6G""ih,f) 



d^^it, s) 



R 



dti 



-G^''{t,ti 



(5*'^(ii) (5*'^(ii) 



ds 



(6.15) 



*(.)=*(-,s) 

Before closing this section we rewrite the adjoint Fokker- Planck operator in the new 
coordinates ^(a)- Prom (6.6) we get easily 



5S 



+ 



G''''{h,t)-^ + a{t) ^ 
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+ 



(6.16) 



It is evident that the additional terms due to our modifications exhaust the most general 
structure of the adjoint Fokker-Planck operator compatible with gauge invariance. 
We finally rewrite upon inserting the special value (6.10) of a as 



dtdti 



5Si 



tot 



+ 
+ 



5^ 



R 



, ^ 



+ Pe{ti,t2) 



(6.17) 



In retrospect the use of the new coordinates ^(ct) has turned out to be successful in order 
to analyze the damping behaviour of the stochastic process. It is explicitly shown that 
the process along the gauge group, although initially giving rise to pure diffusion, can be 
damped by an appropriate tuning of the gauge fixing parameter a. 



7 The equilibrium distribution 

In this section we exploit the remaining freedom for the choice of the positive kernel 7 
and the vector (3^ to define a specific stochastic process of which we are able to derive 
the equilibrium distribution. 

We repeat once more that the legitimation for the freedom of choosing 7 and Pk was 
given in section 5 where we proved that the expectation values of gauge invariant observ- 
ables are independent of 7 and Pk- Therefore, if we are able to find just one stochastic 
process - defined in terms of specific functions a and (3k - where the equilibrium distri- 
bution is known to us, we have succeeded in finding a standard path integral density 
leading to unambiguous expectation values for gauge invariant quantitites. Eventually 
we can compare it with the path integral density derived by conventional techniques, to 
name the Faddeev-Popov procedure. 
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The idea for our procedure is to start from the class of well converging stochastic 
processes (6.15) and choose in a self-consistent way such specific [3^ and 7 that the kernel 
G appears in the Fokker-Planck operators in just factorized form. This corresponds to a 
situation where the fluctuation-dissipation theorem applies, which guarantees that in the 
equilibrium limit the stochastic process becomes independent of G. We will see that the 
equilibrium distribution of our specially choosen stochastic process is proportional to the 
exponent of the gauge fixed total action Stot- 

To become more specific we recall (6.17), where the adjoint Fokker-Planck operator 
was written out in detail. Suppose that we find a solution for (3^ and 7 such that 

G>'''{ti,t2) = I dh[E^k{tuh) + 5^A{.h,h)]5'''[E^t{t2,h) + 8\(5,{h,h)] (7.1) 

then indeed G is appearing just in factorized form and we can simplify drastically the 
adjoint Fokker-Planck operator 



R2 



dtidt2 



S 



G'"'{h,t2) 



as well as the Fokker-Planck operator itself 

, , S 
dtidt2 — 7—^ 
R2 5*'^(ii) 



G^" 



tot 



+ 



(7.2) 



(7.3) 



5*-(t2) 5-^-{t2) 

It is not difficult to give a solution for Pk and 7 which fulfills (7.1). First we deduce from 
(6.13) the conditions 



dhE\{t^,h)5'' E\{t2,h)+/3e{t2,h] 







/i = l,2,3 



as well as 



dU 



E\{h,h) + hitiM) 5'' E\{t2,ts) + Pe{t2,h) =7(^1,^2 



(7.4) 



(7.5) 



and easily find the self-consistent solutions 



Pk{tl,t2) = -E\{ti,t2)+SHe^4t2,h) 



(7.6) 



and 



7(^1,^2)= / dhe^i{h,ti)5ue^i{h,t2). 
Jr 



(7.7) 
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In the last section we claimed that G could be interpreted as a metric. Prom (7.1) and 
from (7.6-7.7) we see that this is indeed the case: We introduce the vielbein E 



with the help of which G is explicitly decomposable as 

G'^"(^i,^2)= / dtsE>^k{ti,h)S'''E''e{t2,ts) 
Jb, 



(7.8) 



(7.9) 



and by construction is a positive matrix. 

We would like to remark, however, that the vielbeins E are anholonomic as in specific 



8% 



7^0 



(7.10) 



so that there does not exist a coordinate sytem with respect to which the vielbeins E 
could be interpreted as the jacobians. 

After this discussion of the self-consistent solutions (7.6-7) for 7 and Pk which lead 
to factorized Fokker-Planck operators we now are interested in writing down the corre- 
sponding Langevin equations. They read 

SStot , SG-'^ithy 



I dh 



-G'"'{tti 



+ E^k{t,h)dW''{h,s)} 



ds 



(7.11) 



5^^(ti) 

*(.)=*(-,s) 

and are of a similar form as discussed at the end of section 3). Like in the case of (6.7) the 
vielbein E factorizes in a gauge invariant part and a gauge dependent rotation matrix 
(which could in principle be absorbed, due to the previously mentioned properties of 
Ito-Langevin equations) . 

We now derive the equilibrium distribution of the stochastic process described by the 
above Langevin equation or by the Fokker-Planck equation with Fokker-Planck operator 
given by (7.3). We remind that we restricted ourselves to the class of well converging 
stochastic processes, see (6.15) as well as (7.6-6.7), so that the Fokker-Planck probability 
distribution is normalizable. Most crucially we have that G is positive, see (7.9); it is 
appearing in the Fokker-Planck operator in factorized form, see (7.3). As a consequence 
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the formal stationary limit of the Fokker- Planck probability distribution can be identified 
with the equilibrium limit and reads 

f^ — Stot 



(7.12) 



Now we use the explicit forms of the total actions to explicitly present our result for the 
path integral density. 

In this paper we discussed two specific gauge fixing surfaces; for the case of the 
variables "^{ai) corresponding to the gauge fixing surface ip^ = 0, after integrating out 
"^f^-^^ and ^('o-j), our path integral density obtains 



(7.13) 



Here the effective action in terms of ^'^ and \E'^ is given by 



dt-- 
2 



1 



(o-i) (o-i) 



1 - 



(o-l) (o-i) 



1 + m 



(7.14) 



and is agreeing nicely with the result of |T3 



For the coordinates ^E'(o-2) the integration over ^(0-2) ^'^^^ leave any trace 

and the path integral measure is just the exponent of the extremely simple effective action 

1 



dt- 
2 



+ + + • (7.15) 



It is evident that this gauge fixing leads to a quadratic action and it is therefore completely 
trivial to calculate the only relevant gauge invariant quantity, the two point function 



lim((v&U)(t))^ + (*L)W)')(^)=2vr 



(7.16) 



It nicely agrees with the corresponding result in Dirac's scheme (see e.g. [IB] and |n 
with the path integral results in |]13| . 



or 
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At the end of this section we transform the final Langevin equation (7.11) back into 
the original cartesian coordinates Amusingly we end up with a Langevin equation 
with a similar structure as in (7.11): First we have with the Ito rules 

d¥{t,s) = [^dti\e\{t,ti)d^'^{ti) 

5eV(t,ti), 



R 



+ 



dU 



Due to the fact that the $ are fiat coordinates we have - similarly as (3.9) 



/ 



dt,^-&^^0 



(5$^(ti) 



so that we finally arrive at 

d¥{t,s) = / dti 

JR 

+ e\{t,h)dW\h,s)} 
Here we defined new vielbeins e 



5$*^(ti) 5$'=(ti) 
$(.)=*(-,s) ■ 



ds 



(7.17) 



(7.18) 



(7.19) 



e'kih, t2) = 5\5{ti -t2)+ dt^e^ih, t^) [-E\{ts, ts) + ^^^6^4(^2, ^3)] (7.20) 

Jr 

and a metric g 

g'\hM)^l dt3e\ih,h)5'"'e'mit2,t3). (7.21) 

"'R 

We recognize that the Langevin equation (7.20) indeed has the same form as (7.11); the 
same conclusions apply about the application of the fluctuation dissipation theorem and 
the independence of the stochastic process of g in the equilibrium limit. 

We close with the remark that although the Langevin equation (7.20) is defined with 
respect to the original cartesian field coordinates the procedure of stochastic gauge fixing 
has introduced a dependence of the new coordinates ^(o-). 



8 Geometrical aspects of the generalized stochastic 
gauge fixing procedure 

The main virtue of the generalized stochastic gauge fixing procedure lies in the fact that 
for specific choices of the extra functions a and (3k the fiuctuation-dissipation theorem 
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can be applied; as a consequence the equilibrium distribution can be derived quite easily. 
Thereby we were led to introduce in (6.13) a positive definite, matrix- valued kernel G'^'^ 
which can be regarded as the inverse of a particular matrix, denoted by Gj^v, on E^^.) x ^. 
What can now be said about its geometrical meaning? 

We will begin with the explicit calculation of G^i,. First it can be shown by a direct 
calculation that the inverse vielbein of E^k which we denote by e^^ is given by 

e'^{t,t') = e\(t,t')+ f dhdhdhZ\tM)K-\h,t2) ■ 

[£^,'(^2,^3) - 5,kZ\h,t2)y^{h,t'). (8.1) 

Here the kernel K^^ G C^(R^, R) is uniquely defined by 

/ dt'K-\ti,t')G^4t',t2) = S{ti - t2), (8.2) 
Jr 

where G^^, denotes the metric inverse to G^", 

G^,(ri,r2)= f dte\(t,n)5ije\(t,T2). (8.3) 
and the following identities hold 

/ dte\{t,,t)E^m{t:t2) - S'^5{t,-t2) 

/ dte%(h,t)E\{t,t2) = (5/5(ti-t2). (8.4) 
The components of the metric G — Jji2 dTidT2Gfj,i,{Ti,T2)S^^{Ti) (8) S^'^{t2) are then given 

by 

Gu.u{ri,r2) = / dte'f,{t,Ti)5ije\{t,r2) 

= / drdr' [ dt^^{t,n)E^\{T,t)GxMr,r') [ dt' E^^ ^{r' ,t')e\{t' ,T2) 
Jb? Jr Jr 

= G^.(ri,T2) + C5^X-i(Ti,T2). (8.5) 



In the above equation we introduced 

G^^in, T2) = G^,{n,T2) - [ dtdt'G^^n, t)K-\t, t')Gi,{t', T2) (8.6) 

Jr"^ 
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Prom a more geometrical viewpoint, the matrix valued kernel J-^dtE^ i{Ti,t)e^ ^{t,T2) 
in (8.5) can be regarded as component of a vector valued one form e on ^[a) ^ Q, namely 

e = j^^ dtidt2dhE\{h, t2)^^{t2, t3)5^^(t3) ® (8-7) 



Insertion of (8.1) yields 

6 



e = J^dt5^^it)^j^- J^^dhdt2dhK-\h,t2)G,^{t2MS^'^^^^ 

+ j^^ dt,dt2K-\t,, t2)5*'(^2) -^^^ (8.8) 

so that for each vector field X on x ^ we finally get the formula 



Here both 



and 



6{X)=X-ZBix) + Zeix). (8.9) 



B{t)=[ dtidt2K-\T,h)G4^{tut2)6^''{t2) (8.10) 



e(r)= / dtiK-\r,ti)5^\ti) (8.11) 

are Lie ^-valued one forms on E^^-) x ^. Since the form B is ^-invariant and B{Z^) — ^ 
V ^ e Lie Q, it defines a connection in the trivial principal ^-bundle E(ct) x ^ — > E(o.). 
However, the corresponding horizontal subbundle, which is given by all those vector fields 
X on S(o-) X Q which are annihilated by B, i.e. B{X) = 0, is orthogonal to the gauge 
orbits with respect to the original metric Gfj,^. This can be proven by a direct calculation. 

Hence the term X — Zb(x) in (8.9) precisely gives the projection of X onto this 
horizontal subbundle. In summary, the global form of G is given by 

G{X, Y) = G{e{X),e{Y)) = G{X - Zb^x), Y - Zb(y)) + G'(Ze(x), ^e(y)) (8.12) 

for any two vector fields X, Y on E^g.^ x Q. From our discussion it is evident that the terms 
appearing in the second part of equation (8.12) correspond exactly to G^;/(ti,T2) and 
5^lA^vAK~^{Tl, T2) of (8.4) respectively. Hence G^^ is the induced metric on the horizontal 
subbundle, which descends to a well-defined metric on E(o-). In fact, for /i, P = 1, 2, 3 this 
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horizontal metric G^p can be shown to be the inverse of the kernel GP'" appearing in the 
Langevin equations (6.7) of the gauge invariant variables. 

Before going on, however, we recall that due to the result (7.10), no coordinate system 
exists such that e could be interpreted as a jacobian. 

The curvature of the connection is given by 



= ^ [ dtidt2 
2 Jr2 



S^^'iti) A5^''{t2). (8.13) 



An explicit calculation yields that it is nonvanishing. Hence there does not exist a mani- 
fold whose tangent bundle would coincide with the horizontal subbundle defined by the 
connection B. 

At the end we would like to give some comments on the geometrical meaning of (7.14) 
and (7.15). The action TZ in (4.2) induces a free action TZ on £ 

n{0,gm = {R{g{t)Mt),^'{t)-g{t)), (8.14) 

in a natural way. As a consequence, the projection n : £ ^ JH := £/Q, with 7t{0) = [<^]a^, 
where ['^]m denotes the equivalence class of with respect to the action TZ, also admits 
the structure of a principal ^-bundle over Ai, where the quotient is taken with respect to 
TZ. Moreover, S Xg A can be regarded as an associated fiber bundle over A4 with typical 
fiber A and projection tt{0, A) = Tf{0) = [0]m- The section ai which was defined in (4.4) 
induces a global section ai of tt : £^ ^ by (Ti([(^];k) = {R{{p^){p,0). Let S(o-i) = im d^i 
denote the corresponding gauge fixing surface in £, then the projection S(o-i) — * ^(o-i) 
yields a fiber bundle with typical fiber A equivalent to n : S Xg A ^ M.. Integrating 
out ^(o-j) corresponds to fiber integration along the gauge group in the principal bundle 
£ X A ^ ^(o-i) and integrating out ^(^-j) corresponds to fiber integration along the fiber 
A in S(cr^). 

The geometrical interpretation of (7.15) is slightly different. The action TZ in (4.2) 



also yields an induced action of G on A, which although it is not free, also ^ gives 
rise to a smooth quotient Af := A/Q. Since every gauge potential A E A can be gauged 
away by a suitable gauge transformation, the orbit space Af consists only of the point [0]. 
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Like in the former case, £ Xg A can be regarded as a fiber bundle over J\f = {[0]} with 
typical fiber S. Integrating out ^^0-2) ^'^^ corresponds to fiber integration along Q in the 
principal bundle £ x A —> S^^a)- Since ^{a^} — {0} ^ ^) integrating out ^(o-j); ^ = 1) 2, 3, 
corresponds purely to an integration over £. 

9 Outlook 

In this paper we performed the stochastic quantization of the helix model and were able 
to show agreement with the path integral procedure in a nonperturbative way. 

We generalized the stochastic gauge fixing procedure by not only introducing an 
extra drift term horizontally to the gauge orbits, but by also modifying the Wiener 
process itself. This modification was done in such a way that expectation values of gauge 
invariant objects remained unchanged. We gave a detailed analysis of the underlying 
geometrical structures of the model. Although the newly introduced extra terms at first 
sight seemed to render our approach complicated, just the contrary was proven to be the 
case: Due to a delicate interplay between the total drift term and the modified diffusion 
term the fluctuation dissipation theorem could be applied and the equilibrium limit of 
the stochastic process was derivable in a straightforward manner. 

As one out of several immediate applications of our nonperturbative stochastic ap- 
proach we are presently studying the helix model for the case where gauge fixing surfaces 
can be defined only locally; another work in progress is the phase space formulation of 
the helix model and the introduction of the BRST scheme. 

Acknowledgement 

We thank H. Grosse, H. Nakazato and K. Okano for valuable discussions. H. H. is grateful 
for the kind hospitality offered to him by the High Energy Physics Group of Tokuyama 
University, Tokuyama and by the High Energy Physics Group of the Dept. of Physics of 
Waseda University, Tokyo, where this paper was completed; he is thankful for financial 



38 



support from the Japanese Society of the Promotion of Science. 



Appendix 



For the case of the variables ^(o-i) (see (4.15) - (4.16)) the vielbeins E and e (see (3.4)) 
are given by, respectively 

^ cosip^{t) -siinp^{t) -*2(i) 0^ 
sin (p^{t) cos (p^{t) -^^t) 
-dt 1 

-10 



S{t - 1') 



I 



( cos (p^{t) sin cp^{t) -ip^it) ^ 
- sin cp^{t) cos cp^{t) ip\t) 



V 













-1 

1 -a 



s(t - 1') 



Furthermore we have 

( 1 -*2(i) q \ ( cosip^{t) -sinip^{t) o\ 

1 *i(i) 

-dt 1 
0-10 



\ 



J 



sin ip^{t) cos (p^{t) 
10 

1 



(A.l) 



(A.2) 



5{t - 1') (A.3) 



In the case of the variables ^(0-2) (see (4.17) - (4.18), (4.6)) we have 

( cos F{t)5{t-t') - sin F{t)S{t-t') -^\t)h{t') -^If\t)E{t,t') \ 

sin F{t)5{t-t') cos F{t)5{t-t') ^\t)h{t') ^\t)E{t,t') 

5{t - 1') - h{t') -E{t,t') 

\ -h(t') -E{t,t') J 



where 



E{t,t') = e{t-t') - / dTh{T). 
Jt' 



(A.4) 
(A.5) 
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Furthermore 

^ cos F{t)S{t-t') sin F{t)S{t-t') 

- sin F(t)5(t-t') cos F(t)5(t-t') ip\t)5{t-t') 

5{t - 1') - h{t)h{t') -S{t-t') 





as well as 



-mt - 1') J 

(A.6) 



( S{t-t') -^^it)hit') -^^it)E{t,t')^ 

5(t-t') ^\t)E{t,t') 

- - -s(t,t') 

J 

( cos F{t) - sin F{t) o\ 

sin F{t) cosF(t) 

10 

1 



(A.7) 
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